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Experimental detection of entanglement of an arbitrary state of a given bipartite system is crucial 
for exploring many areas of quantum information. But such a detection should be made in a device 
independent way if the preparation process of the state is considered to be faithful, in order to avoid 
detection of a separable state as entangled one. The recently developed scheme of detecting bipartite 
entanglement in a measurement device independent way [Phys. Rev. Lett 110 , 060405 (2013)] 
does require information about the state. Here by using Auguisiak et al.’s universal entanglement 
witness scheme for two-qubit states [Phys. Rev. A 77 , 030301 (2008)], we provide a universal 
detection scheme for two-qubit states in a measurement device independent way. We provide a set 
of universal witness operators for detecting NPT-ness(negative under partial transpose) of states in 
a measurement device independent way. We conjecture that no such universal entanglement witness 
exists for PPT(positive under partial transpose) entangled states. We also analyze the robustness of 
some of the experimental schemes of detecting entanglement in a measurement device independent 
way under the influence of noise in the inputs (from the referee) as well as in the measurement 
operator as envisazed in ref. [Phys. Rev. Lett 110 , 060405 (2013)]. 


I. INTRODUCTION: 

Entanglement is shown to be a resource for quantum 
information proccesing, like quantum cryptography [1], 
telepotation [2], quantum metrology [3], channel capac¬ 
ity [4, 5]. Moreover, it helps to speed up computation 
sometimes [6] over the existing classical algorithms. Also, 
understanding entanglement is neccessary to reveal vari¬ 
ous non-classical nature of the physical world. Entangle¬ 
ment in quantum theory correlates two parties in such 
a way that at the individual levels, they lose their inde¬ 
pendent state descriptions. Mathematically, two parties 
A and B are separable iff their joint state p can be ex¬ 
pressed as ® PB^ acting on the Hilbert space 

T~Lab = 0 TLb 1 where aj = 1 and 0 < < 1 

Vj , PA^ and pB^ are density matrices acting on TLj^ and 
Hb respectively. If the joint state can’t be written in the 
above form, we call them entangled [7]. 

Deciding whether a unknown state of a given bipar¬ 
tite quantum system is entangled or separable, remains 
as a challenging problem [8-10] right from the initial 
stage of quantum information theory. People tried to see 
the presence of entanglement in a state experimentally 
through some Bell-inequality violation, but there are en¬ 
tangled states which do not violate any such inequality 
[11]. Given any entangled state PAB, it is (in principle) 
possible to find out an entanglement witness observable 
Wp, whose measurement can distinguish the given entan¬ 
gled states from all possible separable states of the sys¬ 
tem, and the measurement is, in principle, implementable 
experimentally, even using local measurement settings, 
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(possibly) supported by classical communication. It may 
be noted that an entanglement witness is a hermitian 
operator W : Ha ® 'Bb Ha 0 Hb such that there 
exists atleast one entangled state pab on Ha 0 Hb, for 
which TT\W pab] < 0 while Tr[LE ct^b] ^ 0 for all separa¬ 
ble density matrices cjab on Ha ®Hb- So W witnesses 
the entanglement in the state pab- But TT[WaAB] ^ 0 
does not necessarily imply aAB is separable [12]. In gen¬ 
eral, Wp depends on the the form of entangled state pab, 
and hence it does not have a universal character. Thus 
a state-independent witness operator is a desired one for 
most practical purposes. 

The conventional entanglement detection strategies are 
based on the local measurement on each subsystem A and 
B. In this case, perfectness of both the local measuring 
apparatuses is necessary. In fact, if some measurement 
results get lost or some over counting occurs in the mea¬ 
surement, it may lead to an erroneous conclusion about 
the entanglement of the state pab- In ref. [13], time shift 
attack has been used experimentally which causes identi¬ 
fication of a separable state as an entangled one. Such a 
scenario can be avoided if one can witness entanglement 
in a measurement device independent way. 

Recently, Branciard et al. [14] demonstrated im¬ 
plementation of any entanglement witness operator in 
a measurement-device-independent way using Buscemi’s 
work on semi-quantum nonlocal game [15]. This scheme 
of Branciard et al. [14] guaranteed that no separable 
state will be detected as an entangled state, irrespec¬ 
tive of the kind of noise effects present during the mea¬ 
surement [14]. But their scheme depends on the form of 
the shared entangled state. Hence their witness operator 
lacks universal character, that is, it is unable to detect 
entanglement in a unknown state of the bipartite system. 

Augusiak et al. [16] showed the existence of a univer¬ 
sal entanglement witness operator that can detect entan¬ 
glement in any two-qubit state, depending on the Peres- 



2 


Horodecki PPT criteria [7, 17]. But this does not address 
measurement device independent implementation of the 
witness operator, and moreover it is restricted to two- 
qubit case only. 

In this paper, we claim that this two-qubit universal 
entanglement witness operator of Augusiak et al. [16] 
can be implemented in a measurement device indepen¬ 
dent way. We then generalize this idea for higher di¬ 
mensional bipartite NPT entangled states and hnd out 
a set of finite number of universal NPT witnesses for 
any given bipartite system. We conjecture that no uni¬ 
versal entanglement witness operator can exist to detect 
entanglement in a PPT entangled state. We also point 
out that, in the aforesaid scheme of Branciard et al.’s 
measurement device independent entanglement witness 
scheme [14], if some general type of noise is added with 
the inputs or with the measurement operator, the referee 
will not face any problem in witnessing entanglement pro¬ 
vided he knows the character of the noise (here we will 
provide analysis for a particular type of noise and general 
treatment in appendix VII A) . In case the referee doesn’t 
know, we show how robust the witnessing operation can 
be dealt with. 

Our article is organized as follows. In section II, we 
describe measurement device independent (MDI) scheme 
of entanglement witness (EW). Section III is devoted to 
MDI implementation of the two-qubit universal EW of 
Auguisiak et al. [16]. IV describes universal witness 
for NPT-ness of two-qudit states. Section V, describes 
the conjecture about the impossibility of the existence of 
universal witness operator for bipartite PPT-entangled 
states. In section VII A, we analyse the possible noise 
effects in measurement and inputs. We draw the conclu¬ 
sion of our work in section VII. 


II. MEASUREMENT DEVICE 
INDEPENDENT(MDI) SCHEME OF 
ENTANGLEMENT WITNESS (EW) : 

This scheme starts as a cooperative game played be¬ 
tween two parties Alice(M) and Bob(S), who receive 
quantum states {r^ € T~Lao} respec¬ 

tively as inputs from a referee, who then ask them to 
produce outputs a and b respectively. Here 
= dim('H^)= (Ia and dim(HB£,) = dim('HB) = dg and 
for our analysis, we consider (Ia = ds = d. During the 
game, Alice and Bob aren’t allowed to communicate any¬ 
way, but before starting the game they have to decide 
their strategies where they may share apriori a bipartite 
state PAB- If file state pab is entangled state they can 
always get an average pay-off which is more than the case 
when Pab is any separable state. So, from the average 
pay-off values, the referee can conclude whether the cor¬ 
responding shared state was entangled or separable ( by 
knowing the maximum average payoff for all possible sep¬ 
arable states). To produce outcomes of the game, each 
party has to perform local joint projective measurement 


Bell state Bell state 



FIG. 1: (colour online) Alice and Bob performs Bell state 
measurement on the states of the joint systems AoA and 
BBo respectively. Bine arrows are lossless qnantnm chan¬ 
nels for sending input states to Alice and Bob. Red arrows 
correspond to the lossless classical channels for receiving the 
outputs. Green ellipces indicate joint Bell state measurement 
performed by the players. 

on their individual inputs and their individual subsys¬ 
tems of the shared state, see Figure I for details. 

The authors of [14] described the aforesaid projective 
measurement in the maximally entangled state 

- d—1 

l'^}AaA = l'^}BBa = E 1-^^) ^ ® 

with probability of output (a,b) = (1,1) for given 
quantum input pair (Ts,ujt} being Pp^^{l,l\s,t) = 

® I^)bbo (®l)('rs ® PAB ® wt))’ and 
showed that if the average payoff function, Ti{pAB) = 
Y.8,tPs,tPpj,B{^A\s,t) then 

I{pab) '■= max{n(cr_4B)} - n(p_4B) 

= { E (1. l|s> t)| - ^ Ps,tPp^B (1> Ms, t) 

s,t s,t 

= J2Ms,tPp^,{i,Ms,t). 

s,t 

( 1 ) 

where, the maximization is performed over all possible 
separable state ct^b of the bipartite system. The out¬ 
come (1,1) is nothing but the successful projection on 
the state 0 pab G) iot by the tensor product of the op¬ 
erator ($1 <8 \^)bBo and /3s.t is pay-off value 

corresponding to the input state (Ts,uit). 

If we want to use I {pab) as a signature of entangle¬ 
ment witness, it has to be proportional to Pt:\W pab\, and 
the witness operator W needs to be expressed as a lin¬ 
ear combination of the tensor product of the transposed 
input states, 

W = ^Ujf, (2) 

S,t 
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which, in principle, is possible in this case as the inpnts 
Tg (to Alice) and ujt (to Bob)-supplied by the referee- are 
sufficient enough to span the respective vector space of 
witness operators on Ha ® Hb- Then the referee will de¬ 
termine the values of Ps,t accordingly prior to the game, 
such that he can identify whether the shared state is en¬ 
tangled from the sign of I{pab)- Indeed this can be done 
as W is hermitian and transposed input density matri¬ 
ces are also hermitian. The referee just have to choose 
the input states and Wj in order to form a complete 
set of operators { tJ ® wf } to express IT as a linear 
combination of these operators. 

In this scenario, accuracy in preparation of the input 
states must be trusted. In ref. [14], it is proved that no 
noise introduced in the local measurement operator can 
lower the bound I[pab) ^ 0 for separable pab^ so, in 
this scheme, separable state will never be identified as an 
entangled one. In this sense the scheme of ref. [14] is 
measurement device independent. 


A. Detecting entanglement in unknown 
two-qubit state : 

In the case of witnessing in a unknown two-qubit state 
[16], instead of sharing single copy of the state, the play¬ 
ers need to share four identical copies of the state. This is 
so because the signature of det(p^g) = product of eigen¬ 
values of = Tr[lTp®g] = A 1 A 2 A 3 A 4 , completely de¬ 
termines whether the two-qubit state pab is entangled or 
separable. Here A^’s are eigenvalues of p^g and IT is a 
hermitian operator acting on H^ 0 H^^ (with dimHA 
= dim'HB = 2). From the Newton-Girard formula [18], 
we have : A 1 A 2 A 3 A 4 = 

l-6^Af-k8^Af + 3(^^A^^ - 6 Aj ^ 

i—1 i—1 i—1 ^ 



To formulate this in terms of an operator, we first use 
the swap operator T^^i, defined as 


= l‘^fe)AiBi ® I‘/’i)a2B2 ® (3) 

with the property Tr(T('=)p®'') = Tr(p^) [16, 19, 20] for 
any hermitian matrix p and T^^^ is not, in general, a 
hermitian operator except for k = 2. So to make it an 
observable in place of T^^^ 5 ( 1 ^^^^ + T^^^I) will be used 
for farther analysis. As the swap operator is a real 
matrix with respect to a fixed basis, used for the whole 
analysis. 



is given by : 


= Tr(TW(p5-g)®'=) 

= iTr((T('=^ -k 


y{k) 0 y(k)T y{,k)T ^ yik) 

where T^'^^ : HJ'" 


«). 


® \4>2)j^ 0 .® \4>k)j ^) 

= \(t>k)j^ 0 0 .0 \(l)(k-i))j^ for j = A, B. (4) 


Then, hermitian operator witnessing entanglement in 
arbitrary two-qubit state, is given by. 


J 


q/univ ^ _ i (^A) 0 vA)T + yA)T ^ 

+ ly(2) 0 V(2) 
8 

and our I [pab) turns out to be proportional to 
Tr(IT'^“^p®g) = det( p^^), see the next section III. 


The above form of witness operator doesn’t depend on 
the type of two-qubit shared state, so that entanglement 
in an unknown state can be detected. In this sense it’s 
universal entanglement witness operator. 


+ ^I4x4 ® ® + T(3)T 0 y(3)^ 

— ^IIl6xl6 G) T^^\ (5) 


III. MDI IMPLEMENTATION OF TWO-QUBIT 
UNIVERSAL EW : 

Recent work by Bartkiewicz et al. in [21] demonstrated 
the implementation of the aforesaid universal entangle¬ 
ment witness operator (5) for two-qubit case, using pho¬ 
ton polarized states. But that has not been done in a 
measurement device independent way. 

Authors of ref. [22] and [23] , address entanglement de- 
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tection of arbitrary two-qubit states in MDI way, but 
their approach focus on optimization of entanglement 
witness operator instead of the question of the univer¬ 
sality of that, as we discuss throughout this article. 

The game of Fig 1, starts with two inputs (one for Al¬ 
ice and one for Bob) and single copy of a shared state, 
while our universal measurement scheme for witnessing 
entanglement needs four identical copies of shared state 
PAB S't each run of the experiment and the witness op¬ 
erator is a 256x256 dimensional matrix. Thus, 

the input Hilbert space is ^ [c 2]®4 ^ ^ 

instead of ® C^. As a result of that, both Tg and w* 
must be density matrices on 

So, the witness operator of equ. (5) should be 

of the form. 


= ( 6 ) 

S,t 

where, the input states and uJt are expandable in ap¬ 
propriate Gell-mann matrix basis such that 


A.Tr(rjA) 
® 16 

rp A.Tr(a;^A) 


= Ea. 




Wt = 


16 


-Ea, 


-3 9 


(t) 


(7) 


where A is the vector of the 16 x 16 generalized 
Gell-Mann matrices [24], whose first component repre¬ 
sents identity matrix. Here, i,j € {0,1, 2,3,..., 255} 
and Ai, Aj are i,j—th components of A = 
(I, Ai, A 2 , A 3 , ...,A 254 , A 255 ) respectively. Also 
and are real coefficients. 

Then from the above relation ( 6 ) and using the actual 
form of the universal witness operator (5), the referee will 
calculate the quantities. 


Tr(lF""-(A, 0 A,)) = (16)2 qf\ 

S,t 


that will give a set of (256)^ number of linear equations 
for all i,j, from which the referee calculates the coeffi¬ 
cients /3s_t in equ. ( 6 ). Also, instead of using a two- 
qubit Bell state projector |(/)+) ((/)+1 = ^ j=o N) (jj l 
or it’s local unitary equivalent as a measurement oper¬ 
ator, each party should use Bell state projector of two 
sixteen dimensional systems as the measurement oper¬ 
ator |$+) (4>+| = I**) 0.11 it’® unitary 

equivalent, for each party. Then the negativity of 


I{pab) = det(p 5 B) = Tr[VF““''(p^iBi 0 PA2B2 ® PA3B3 ® ^^464)] 
(4>+|0| 

^’'’)bo:BiB 2B3B4 (d>+|)( 


0 pAiBi G) PA2B2 ® PA3B3 ® PAiBi ® J 


( 8 ) 


confirms the presence of entanglement of the shared state. 

Here if we choose to work in the computational ba¬ 
sis {jo) , 11)1, any pure state of a single party A, can be 
written as, 

\M = {a |0)i -h Pi |1) 0 (02 |0) -f P 2 |1) 

0(03 |0) -h/33 |1))^3 ® (04 \ 0 )+ Pa 
such that the set {ji/’^t), |^b)} spans the support of 

„®4 

Pab- 

The form of the swap operators has to be found out by 
their action. 


Therefore, 


y(4) ^ E (bj) actson'H^®'‘ ( 10 ) 


1/(3) 0 I= (ag |0) + /33 |1) 0 {ai |0) + /3i |1) 

0(02 |0) -h P 2 |1) ® (a4 |0) -h Pa |1) 

( 11 ) 

Therefore, 


yW IV’^) = (aA |0) + Pa |1) 0 (ai |0) + Pi |1) )^^ 

0(02 | 0 ) -h P 2 |1) ® (as | 0 ) -f p^ | 1 ) 

( 9 ) 


y(3) ^ E acts on ( 12 ) 

Similar expressions are valid for the operaors act on 
'Hb®'^ where, n = 3,4. 
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A general support vector in 'Hab^'^ can be written as, 

I'^Pab) = (ai |0) + /3i |1) 0 [a'l |0) + j3{ |1) 

0(^2 |0) + /32 |1) 0 {a '2 |0) + /32 |1) )g2 

0(a3 | 0 ) + /33 ® (ttg | 0 ) +/33 |l))g^ 

0(a4 |0) + /34 |1) ® K |0) + P'i 11 ) )b, (13) 

Therefore, 

^ y( 2 ) 

= (a 2 | 0 ) + f32 | 1 ) 0 (a 2 | 0 ) + | 1 ) )b 2 

0 (ai | 0 ) + /3i | 1 ) 0 (a'l | 0 ) + /3^ | 1 ) )g^ 

0(04 |0) + P4 |1) )_4^ 0 {ot'i |0) + ^4 |1) ) b ^ 

0 ( a 3 | 0 ) + /33 | 1 ) )^3 ® (4 | 0 ) + 4 | 1 ) ) b 3 ( 14 ) 

Thus, 

1 

^ \kl,ij) {ij,kl\ acts on TLab^^ (15) 

IV. UNIVERSAL WITNESS FOR NPT-NESS 
OF TWO-QUDIT STATES: 

Any two-qudit density matrix which has at least one 
negative eigenvalue after partial transposition, is called 
an NPT state. Any bipartite state having NPT is nec- 
cessarily entangled. Now, to determine NPT-ness of a 
unknown two-qudit state (i.e. to determine whether the 
state has NPT or PPT), we need to check negativity of 
the eigenvalues of [25] . Let us look at the character¬ 
istic equation for the matrix p^g : 

flet(4B - >'h^xd^) = 0 , 

which is in the form, ^(— 1 )^^ ~^ad 2 _iX^ = 0 , 

i=0 

( 16 ) 

where, the coefficients oq, ai, 02 ,...., 0^2 can be given in 
terms of the eigenvalues Ai, A 2 , A 3 ,..., A £;2 of pab as fol¬ 
lows : 


d^ 

ao = 1 , ai = Ai, 02 = A^A^, 

i=l 


0.3 


E 


A j Aj Afc, 




0(42 — Aj 
i=l 


As p^g is hermitian, A’s are all real, and thereby Oi’s 
are all real. If Oi > 0 V i, we have pjg > 0. Oth¬ 
erwise Pab has NPT. Now, we can determine whether 
negative roots of the characteristic equation exist. We 


can write each of these coefficients as a polynomial in 
As’(eigenvalues of p^g) using the Newton-Girard formu¬ 
lae [18]. In fact we have : oi = 1 because partial trans¬ 
position operation is trace preserving; 


1 

= 0(1 - E-^^) = 41r(IV2p5|), ’^ilh. 




W2 — - 


(17) 


“3 = ^ E ^ E ^0 = Tr(IT3p®|), 

2=1 2=1 

with, W 3 = - 3 IId 2 xd 2 0 

-b 1/(3) 0 y{3)T y{3)T ^ -j;>(3)j . 


(18) 


and, 04 = 

2=1 2=1 2=1 2=1 
= Tr(W4p®"g) 

(19) 

where W 4 has the same expression as that of in 

Equ. (5); so on and so forth. 

Our argument starts from hnding the signature of the 
coefficients o^ with i > 2. For most favourable case, 
02 may turn out to be negative and thereby pab has 
NPT. So, two copies of the shared state, pab is enough 
to determine NPT-ness while both the inputs Ts and 
udt are density matrices on Hilbert space (C'’*)®^. But 
for the general case, the players have to share b’ num¬ 
ber of copies of the state pab whenever referee wants to 
know the sign of the coefficient o^. Also the dimension 
of the input states (ts and ujt to be supplied by the ref¬ 
eree) should increase accordingly. For the worst case, we 
have to go upto the coefficient 0^2 for which we need cP 
copies of Pab and Tg and ujt will be density matrices on 
the Hilbert space As, for any PPT state paBi 

all of Ai, A 2 , A 3 ,..., Ad 2 are non-negative, all the coeffi¬ 
cients 02 , 03 ,..., a £;2 are non-negative. So, the operators, 
H/ 2 , W 3 ,..., ITd 2 indeed work as witnessing NPT-ness of 
the state pab- 

Note that, these witness operators don’t depend on 
the shared state pab i and so, we can use the scheme 
similar to that described in section HI to implement the 
measurement of these operator in measurement device 
independent way. A notable point is that, for two-qubit 
case, under partial transposition, the shared density ma¬ 
trix can have only one negative eigenvalue, if the state is 
entangled [16]. Thus for the case of two-qubits, the sign 
of the coefficient 04 is enough to decide the NPT-ness, 
and thereby entanglement of the state as done in section 
HI. We don’t need to check the signs of the coefficients 
02 , 03 . But for higher dimensions, that is not the case. 
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A. Comparison between our MDI 
scheme and the conventional state 
tomography : 


In state-tomography, for a two-qudit system, we surely 
need, —1 number of measurement settings, as the num¬ 
ber of coefficients in any basis representation ( like gen¬ 
eralized Gell-Mann matrices A„ [26] ) are — 1 : 


1 

Pd^xd^ = 


d*-l 

I 0 I + ^ 

n=l 


( 20 ) 


Coefficients r„ £ M, can be obtained from the expecta¬ 
tion value of observable A„ : = Tr[ . There¬ 

fore, for each coefficient, atleast single copy of the state 
has to be used. In order to reconstruct the state, we 
need to know exact values of the coefficients, their signs 
only don’t give the sufficient information. The erroneous 
numerical values of these coefficients lead to a wrongly 
estimated state, So, it is not easy to think of a 

scheme like mesurement device independent tomography 
for a two-qudit system. 

But in our case, without going into any MDI scheme, 
to measure 02 in Equ. (17), we have to measure the 
expectation value of a single observable only. For 03 
in Equ. (18), we need to measure, expectation values of 
two observables 1^2 ,,,^2 0E( 2 ) ^ ^( 3 ) ^y(3 )t_^ y(3)T^y(3) ^ 
Similarly for 04 , we have four observables to measure. 
So, we need very less number of measurement settings 
compared to the tomography, at the cost of n—copy usage 
for coefficient a„. 

In the universal MDI scheme of finding NPT-ness of 
any two-qudit state, for the coefficient a„, we need two 
general Bell state projectors in Hilbert space 0 , 

and the input density matrices Ts, ujt on Hilbert space 

2 

C” . These inputs will span the vector space of HAi, so, 
their required number will depend on the number of ob¬ 
servables in the expression of Wm which is in general less 
than the observable in conventional tomography. For our 
case, the sign of the coefficients a„’s are enough to decide 
whether the shared state is entangled or separable. 


V. WITNESSING AN ARBITRARY PPT 
ENTANGLED STATE OF A GIVEN 
BIPARTITE SYSTEM: 

Existence of PPT entangled state (so called 

PPT bound entangled states) is well-known whenever 
dim('H^) X dim('Hg) > 6. Evidently our scheme of uni¬ 
versal NPT-ness witness, as described in section IV, does 
not work here to detect entanglement of any such state. 
Below we try to argue that a universal entanglement wit¬ 
ness operator can not exist to single out entanglement in 
any PPT state of a given bipartite sytem, through some 
conjectures. 

Gonjecture I : If dim('H_ 4 ) x dim('Hg) > 6, there can 
not exist one (or a finitely many) universal entanglement 



FIG. 2: (colour online) Distribution of density matrix spaces 
on Hilbert space Ha (SHb- [ (I) -b (H) + (III) -b (IV)] = 
T>ab, the set of all density matrices. (I) = Hab — Vab is 
the set of all NPT states, Vab is the set of all states having 
PPT (II) = Vab is the set of all PPT entangled states, [ 
(III) -b (IV) ] = Sab is the set of all separable states, [ (II) 
+ (HI) ] = T^'ab i® ihe set of convex combination of PPT 
entangled states. Region (III) clearly depicts the non-empty- 
ness of iVABO^AB)- 


witness operator detecting entanglement in an arbitrary 
PPT state p^g, and which can be realized in a MDI way. 

Reason behind this conjecture is another conjecture, 
Conjecture II : There can not exist a universal EW 
(or, a finitely many EW operators) for all PPT entangled 
states of any given bipartite system 

Let us try to provide some geometrical argument which 
potentially may give rise to Conjecture H. 

Let Vab , Sab j Vab are respectively be the set of den¬ 
sity matrices, separable density matrices, density matri¬ 
ces which are positive semidefinite under partial trans¬ 
position (PPT), on Hab', Vab be the set of all entangled 
density matrices (PPTE) in Vab- Sab and Vab are both 
convex sets. It is obvious that. Sab C Vab C Vab 
and Vab ~ Sab = Vab- Note that Vab is not con¬ 
vex. For max{dzTO('H^), c?im(7^g)} = 3 together with 
mm{dim{'HA) , dim{'HB)} = 2 , Vab = Sab- Let V'ab 
is the set of states formed by the convex combination of 
states G Vab- 

Edge state : pab G Vab is an edge state iff the state 
{PAB + P0-^b)/(1 + p) G Sab for all p € (0,1] and all 
uab G Sab- Edge states lie near the boundary of Sab 
and exist for all dim('H^,'Hg) > 3. 

It is well known that [27, 28] there always exist atleast 
two edge states Pab^Pab ^ Vab such that, their con¬ 
vex combination, pPab + (f “ p)Pab = G Sab 

for 0 < p < 1. Hence for, dim{'HA),dim{'HB) > 3, 
V'ab G Sab ^ ‘k (null set). This is the case for single 
copy. Fig.2 depicts existence of these edge states in re¬ 
gion (H) near the boundary of (H) and (HI). 

From Fig.2 it is clear that, there can not exist a her- 

mitian operator W : 0 —>■ C'’* 0 C"’* such that 
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Tr[VFp^B] < 0 for all £ Vab while Tr[Vl^(T_ 4 B] > 0 
for all cr_ 4 B £ 5 _ 4 b . So, in the single copy case no such uni¬ 
versal witness operator W exists to witness entanglement 
in every p^g £ Vab ■ 

Now, the question is — As in the case of qubits, is it 
possible to have a ( for finite number of) universal witness 
operator(s) which can detect entanglement in any 

PPT state pAB of two-qudits in case n—copies of the state 
are supplied with n > 2 7 

We conjecture that, such a does not exist. 

Let’s consider, 27^„g„ is the set of all density matrices 
in is the set formed by the convex 

combination of the states cr®g, where, cr^g £ Sab and 
~ (”) 

V'ab the set formed by the convex combination of 
the states p®g, where, pab & 'Pab- If Sa„b„ be the set 
of all separable states in Ha^^ ® it is evident 

that, 5^g^"^ £ Sa„b„ and it is most likely that, Sa„b„ LI 

V'ab 7 ^ (null set), for all n with dim{'HA), dimiHB) > 
3, because of the existence of edge states in 'Ha®"' ® 
Hb®". 

Conjecture III: S'ab^"'^ H 'V'ab 7 ^ ^ (null set) for all 
n with dim{'HA),dim{'HB) > 3. 

Fig.2 depicts Conjecture III is true for n = 1. 

That means, there always exist atleast two states 
Pab’Pab ^ Pab such that, the convex combination of 
their n number of copies, 

Pip'Af" + (1 - ( 21 ) 

r 

is a density matrix describes a separable state £ 'Ha®" ® 
Hb®", cr^B ^ ^AB, and this holds for any n £ N. 

Where N is the set of all natural numbers, 0 < Pr < 1 
for all r and ~ I- 



FIG. 3: (colour online) Distribution of density matrix spaces 
on Hilbert space <»'Hb®"■ [ (I) U (H) U (HI) U (IV) U 
(V) ] = T>A„Bni the set of all density matrices. (11)= SAnBn 
is the set of all separable states, (III) = S'ab^"^ is the set of 
convex combinations of all states like crjg, where aAB € Sab- 

(IV) = 'P'^ab is the convex combination of the states like p®)), 
where pab € Vab- (V) contains convex combination of all the 
states like p®g, where pab € 'Vab—'Pab, set of all NPT states 
at the single copy level. We conjecture about the non-empty- 
ness of the intersection (III) H (IV), and the empty-ness of, 
(III) fj (V) = <(> (null set). 


for a complete set of linearly independent observables 
O : 'Ha®" ® Pb®" ^ Ha®" ® Pb®"- Even verification 
of this one may turn out to be difficult. 

This difficulty seems to arise from the well-known fact 
that the separability problem (that is to find out whether 
arbitrary state of a given bipartite quantum system is 
entangled or separable ) is NP-hard [29, 30] . 


Validity of conjecture III implies nonexistence of any VI. NOISE ANALYSIS FOR REAL 

hermitian operator, EXPERIMENTS : 


: Pa®" ® Pb®" ^ Pa®" ® Pb®", 

for which, Tr[lV 4 „g^ (p^g)®”] < 0 for some pab S 'Pab 

together with Tr[lF^„g„ (cr^tg)'^”] > 0 for all cjab S Sab 
— irrespective of the choice of n. 

Thus, validity of conjecture III implies that there can 
not exist a universal EW (or, a finitely many EWs) which 
can detect entanglement in all the PPT (bound) entan¬ 
gled states of ^ -I- S whenever d > 3. 

This atomatically implies the validity of the conjecture 
II. It appears to be quite difficult to verify conjecture 
III directly. One may try to verify whether for some 
given states Pab € Pab, there exist a separable state, 
aAB G 'Ha®" O 'Hb®" such that. 


Tr 


^ J^P^iP-As) 


(a)\®ra 


= Tr 


O aAB 


Recently, a few experimental works [13, 31] have been 
done to implement this MDI protocol to test the entan¬ 
glement of two particular type of two-qubit states. 

We now analyse the effect of possible noises in the Bell- 
state-measurement (BSM) part of their experimental se¬ 
tups. Here we introduce three kinds of possible noise (see 
Fig. 4) : 

(1) Photon loss in the polarizing beam splitter 
(PBS) : We model here the photon loss in a PBS through 
the action of a white noise - as described below. 

(2) Error in the angle of rotation in the half 
wave plate (HWP) : HWP rotates the polarization 
axis of the light vector (which carries the information 
about the projected state by the measurement operator) 
with respect to it’s direction of the propagation of light 
ray by the angle equals two times of the angle (9) between 
the fast axis of HWP and the polarization axis. So, if the 








FIG. 4: (colour online) Bell-state projection ^{\HH) + 

VV){{HH\ + measurement via coincidence counts. 

Photon propagations are denoted by red-rays. 



FIG. 5: (colour online) Light ray passing through Half wave 
plate. The polarization axis, rotates by angle 29 with respect 
to the direction of propagation. 9 is the angle between the 
fast axis of the direction of polarization of the incident light 
ray. The first axis of the HWP lies in the a;z—plane. 


angle of rotation be improper, we can it as 29 + ry, where 
rj is the error in BSM. 

(3) Detection inefficiency in the PBS and Diode 
photon-detector parts : For detection inefficiency, 
the probability of detection sets reduced by a factor ^ 
such that 0 < ^ < 1 and the probability of detection 
become ^P(l, l\Ts,uit). 

For our discussion we have chosen |$+)(<i)+| a Bell 
state projector as measurement operator ( as the pro¬ 
jector of BSM), where |<i>+) = -^{\HH) -|- \VV)) with 
the consideration that the horizontally polarized state is 
\H) = (1 0)^ = |0) and vertically polarized state is 

|I/) = (0 1 )^^| 1 ). 

Under additive white noise (a special case), our actual 
measurement operator will be, yx where 

/i is the corresponding visibility. 


We assume the direction of propagation of photon is 
along y axis (denoted by red rays in Fig. 4 and Fig. 
5), where the reference coordinate frames fixed at every 
HWP. Thus if the angle of rotation for the polarization 
axis (shown in Fig. 5) of one player and that coming from 
the referee are respectively 2gi and 2 g 2 then, rotated 

|$+) = (g) |$+) 

= cos (52 - gi) + sin(g 2 - gi) |^'~) ( 22 ) 

where, I'll”) = -^{\HV)— \VH)). So, if gi and g 2 are the 
same or ± 7 r or ±27r then there is no error in the HWP. 
Therefore after passing through noisy PBS and HWP, 
our noisy Bell state projector will be, 

P = gcos^{g 2 - gi) |$+) ($+| 

-bixcos(g2 -5i)sin(g2 - 5i)( 1^'”) 

+gsm'^{g 2 - gi) (T'”| + 

(23) 

Let’s take g 2 — gi = A. Then the noisy projection oper¬ 
ator shared between A and a part come from the referee 

Ao = 

PaoA = Ml cos 2 (Ai) |$+) ($+1 
-b/xi cos(Ai)sin(Ai)(^ ($+| -b |$+) (^'“| ) 

+Misin2(Ai)|vI/-)(T--| + l^I (24) 
Similarly for B and Bq, 

PbBo = M 2 C 0 S^(A 2 ) |4>+) (4>+| 

+g.2 cos(A 2 ) sin(A 2 )^ ($+| + |<i)+) (T'“| ^ 

+M2sin2(A2)|^-)(4'-| + i^I (25) 

Here, gi and yi 2 are the visibilities of the BBS’s of 
the party A— side and H—side respectively in the noisy 
measurement of |$+). Ai and A 2 are the corresponding 
errors in rotational angle in HWP in A— side and H—side 
respectively. 

So, our desired (modified) quantity will be, 


ImodiPAB) — 

Tr {PAoA® PbBo){'^s® PAB^dJt) 


S,t 




S/A \ 2 /\ \ I Ml(l M2) 2 


P 1 P 2 cos (Ai) cos (A 2 ) -b 


cos"^(Ai) 


. AlnA c„»»(A,) + a-w)(i-M2) 

4 y 

I{pab) + additional term. 

(26) 


The additional term is in Equ. (42) and (43) of appendix 
VHB. 
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Now, for simplicity, let us assume that fJ.i = fJ .2 = fJ- 
and Ai = A2 = A. 


In the case of a two-qubit shared state pab, the input 
states Tg and uJt takes the forms. 


J 


I -I- fg.CT 


Wt = 


l+TfS 


2 ■’ = 4 


CS J A 

^I(g)I-l-^ai I(g) (Ti-I-^ cri®I+^Cij ai®(Tj^ (27) 


i=l 

p 3 


i—1 i,j—l 

3 


with, = {xs,t, ys,t, Zs,t) e R , < 1, (01,02,03) G M"*, (61,62,^3) e 

and (cij)f_j_i is a real 3x3 matrix. 


The corresponding form of I-mod{pAB) is given in, ap¬ 
pendix [VIIB] 

According to ref. [31], the input states are given in 
table -I, 


Input 

states 

Ts,U!t 

Bloch vec¬ 
tors Ts.rt 

\H) {H\ 

( 0, 0, 1) 

\v){v\ 

( 0, 0,-1) 

\D) {D\ 

( 1, 0, 0) 

\D){D\ 

(-1, 0, 0) 

\L) {L\ 

( 0, 1, 0) 

\R) {R\ 

( 0,-1, 0) 


Table - I : Input states used in the ref. [31] 
for Alice and Bob’s side and their corresponding bloch 
vectors. 

The corresponding payoffs are : 

Ph.h = Pv,v = I3 d,d = Pd,d = Pl,l = Pr,r = 

Ph,V = I3v,H = Pd,D = Pd,D = (^L,R = Pr,L = “ c - 

(28) 

Other /3s,t are equal to zero. 

The shared state is, pab = P I't'”) -I- I which 

is separable iff 0 < p < i. Thus in Equ. (27), cn = 
C22 = C33 = -p. Cki = 0; when k ^ = 0 

yk,l G {1,2,3}. 

Using Equ. (28) and table-I we have, 

= 1, 

S,t 

^ ^ Ps,t ^ ^ l^s,t ysVt — ^ ^ Ps,t — I5 

s,t s,t s,t 

Ps^txsvt = y] Ps,txszt = y] ps^tvsxt 

s,t s,t 

~ ^ ^ l^s,tysZt — ^ ^ ^s,tZs^t ~ ^ ^ Ps,tZsyt ~ 0 - (^ 9 ) 

s,t s,t s,t 


Therefore, according to Equ. (26) ( see Equ. (44) in 

Appendix VIIB ), 

3^mod(p.4B) = I [1 “ ~ cos(4A)] (30) 

For the allowed values of the noise parameters p, A we 
have, 0 < Imod{pAB) < Hence for these kinds of errors 
no separable state will never be detected as entangled 
state. 

Ref. [13] deals with the entangled state, p = (1 — 
r) l^--) («'-] -f §(li7i7) {HH\ + \VV) {VV\) with 0 < r < 
1. pab is separable for | < r < 1. Here, for this case, 
cii = C22 = (?' — 1) and C33 = 27- — 1. Cki = 0; when k ^ 1] 
ak=bk=0, VA:,Zg {1,2,3}. 

Here we are going to use the data of the Table H of 
the Supplementary section of ref. [13]. 

for that case. 


s,t 

E f3s,tXsXt = E Ps,tysyt = 1) E = 1, 

s,t s,t s,t 

^ ^ ^ ^ Ps,tZsyt — ^ ^ l^s,t^syt — 

s,t s,t s,t 

^ ^ I^S,t^sZt ^ ^ ^s,tys^t ^ ^ I^S,tysZt 0; 
s,t s,t s,t 

(31) 

Therefore, according to Equ. (26), ( see Equ. (44) in 

Appendix VHB ), 

ImodipAB) = |[(3r - 2)p^cos(4A) -|- p^(r - 1) -I-1]. 

(32) 

For ^ < r < 1 and for allowed values of noise parame¬ 
ters ^,/i. A, we have 0 < ImodiPAB) A So, as in the 
case of Ref. [31], separable states will never detected as 
entangled. 
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VII. CONCLUSION: 

Using the prescription of Auguisiak et al. [16] for wit¬ 
nessing entanglement in a unknown state of two-qubits, 
we have extended here the measurement device indepen¬ 
dent entanglement witness scheme of Branciard et al. [14] 
to a universal measurement device independent entangle¬ 
ment witness scheme for two-qubit states, the caveat be¬ 
ing that we need four copies of the state at a time and the 
referee should supply the input statesfrom sixteen dimen¬ 
sion Hilbert spaces of Alice and Bob. In an aim to extend 
this result to higher dimension, we provided here a mea¬ 
surement device independent scheme for universally wit¬ 
nessing NPT-ness of a unknown state of any given bipar¬ 


tite system - provided many copies of the state is being 
supplied. Of course, that number of copies is much less 
than that required for state tomography. We conjectured 
that, there doesn’t exist a single or a finite number of uni¬ 
versal witness operators for witnessing entanglement in 
an arbitrary PPT state of a given bipartite system. Our 
noise analysis of the Bell state measurement scenario in 
both the experimental demonstration of MDIEW [13, 31] 
are in conformity with the MDI of EW of Branciard et 
al. [14]. 

ACKNOWLEDGEMENT: The authors would like to 
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I 

APPENDIX 


A. Inputs and measurement operators are noisy 

An usual apparatus is generally influenced by ambient noise. In [14], the players and referee trust the preparation 
of the input states. So, the noise in Ts and ojt (if any), is supposed to be known by the referee (or, even by the 
players). Some noise may affect the actual shared state but as we are going to detect it, we will consider the 
noise induced shared state as the actual state pab, need to be detected. 

Any general d—dimensional noise can be expressed as ^^.Adxd where n* and A® are the components of generalized 
Bloch vector and Gell-Mann matrix in d-dimension respectively, A° = Idxd [26]. If the visibility p, of the actual input 
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state is same for all input states, and assume the visibility remains constant throughout the experiment, then the 
players (Alice and Bob) would receive the input states. 


Ts = (J-Ts 


(1 -M) 


rii ■A.dxd — 


(1 - m) ^ , r V 

—-—m + -ly[TsAdxd) 


■A. 




and 


(1 - /i) ^ , fJ-rj. f y V 

— 2 —^2 + -ly(uJtAdxd) 


■Ad 


Xd 


instead of Tg and ujt respectively . If the referee does not aware about the noise mixed with the inputs, there will be 
mismatch of the value of I{pab) from it’s actual value if the refree knows the noise. 

In the game, the form of W is known to the referee, he calculates the pay-off functions, Ps,t according to the inputs 
[2]. If referee uses Tg and uit as the inputs instead of t' and his calculated Ps,t will be different than the actual 
that will give the actual value of I{pab)- 
Modified witness operator looks like. 




tT, 


'T , 


'T 

UJ, = 


<F 


s,t 


T T 

:,tTs - 


additional noise terms. 


(33) 


So, by knowing the amount of noise induced in inputs, referee can calculate the modified value of /3g_t, which comes 
from the linear equation generated by, Tv\W'{t^ 0 for all s,t. 

So, we can claim that if the referee knows about the character of the noise in inputs, entanglement determination 
shouldn’t be erroneous. 


In general, the measurement operators can also noisy, and the character of the noise is not supposed to be known 
by the referee nor by the players. But that noise can obviously affect the final decision based on the value of I{pab)- 
If a general additive noise, ^ is added with the maximally entangled projector i I**) Oil > then with 

the visibility v of the original measurement operator, the actual measurement operator appears as 

i' A — y) ^ v 

^ m OjI h—^2— 

i,j 

provided this is a positive operator, that is V |x) we have. 


^^1 E 1*0 Oil + ^ 2 ’^^ rn.Ad^^d^'] lx) > 0 , 

^ ^,3 ' 

In the standard basis,{|i, j)}, we can write 

^ {a,b\A\p,q)\a,b){p,q\ 

a,b.,p,q 


For notational convenience, in place of PABy we will be using p here and also A for A^ 2-^^2 . 


Pp(l,l|s,t) = Tr 


(3 E 1*0 Oil + ® E 1****^ ^*'*^1 + ^^^^2^*'*2.A 

'"IT ' -j j 11 


Ts® P'iPUJt 


V 


(rj p 


i.(l - y) 


(P 


Tr 


(^Tg^ 0 Ttbo [(to 2 .A)(I 0 Wt)]^ p 




(^Tr^o[(Tg 0l)(?7ii.i^)] 0 wf) p 


( 1-^0 

p 


Tr 


^Tr^o [(Tg 0 1)7711.A] 0 Ti'Bo [ 7712 .A(I 0 wt)]^ p 






( 34 ) 


If in (34) toi.A = 7112 .A = I^^d 0 Id(gid , will be denoted here as I 0 I, and if define A = (I, T), then 
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U) 


rp ^ 1 ^ p 


1 /( 1 -z/) 

H - tVAo 


{Ts 0 I)(l0 I) 


0 ujf ■ 


{l-vf 


<p 


Ts ®trBc 




tTAc 


(Ts 


Tree 




) LU. 


p z/(l — Z/) 


1^(1-u)^ 


I + ^ .. 0 ujT ■ 


(I 0 l)(l 0 a;t) 

I (g) I(I 0 Wt) 

{1-^f 


d3 « ^ - t ' ^4 

,2 


I0 I 


Jjl 0 I + ^ (l 0 f) ^ ® ^ (fir(fTj)) 0 (fir(fa;t^)) 


= ^ (I + .tr(frj)) (I + vf.trifujf)) 


(35) 


Thus this is just a shrinking operation of the generalized Bloch vector T. 

If the referee is unaware about the noise then, there is deviation from the actual results.So, referee will use the /3 
functions same as that obtained without noise.In presence of noise in inputs and measurement operator,the witness 
operator form becomes. 


+ 


d4 


trAc 




z/ 


^J-Ts 


d) 


tTi.A 


dxd^ ® (, 


fiUJt 


(1-m) 


ri^.A 


V 

dxd j 


Ki-^) 

iy{l - zz) 


( 1 -m)-, r 

fJ.Ts-\ - - — ni.Adxd] ®trBo 


) ®trBo (m2.A)(l0 (/rwz + E_E.^2.Adxd)) 


d3 


tTAc 


((m'Ts + -^l-^-^ni-Adxd) ® I)(mi.A) 0 (/rwz + ^ ^.Adxd^ 


d 

(1 - d) 


(1 — /r) r 

M'Ts H-1— ni.Adxd 


:) 01 ) 


mi. A 


) trg. 


m 2 .A (10 ( 


(1 — /r) ^ 7 * ,. 

/rwt H- - — n2.Adxd j 


= W""(say) 
(36) 


Therefore, to get the correct result, in the expression of IT, given by the equation [2], instead of tJ 0 ujJ referee 
should use 




Kl-^) 

d3 


Tr 


Ao 


(ts 0 I)(mi.A) 


■ 


d^ 


) UJf 


TiUc 


iz(l — v) rp 

d3 

(ts 0 I)mi.A 


Trso 
0 Trgo 


(m2.A)(l0 wt) 

m 2 .A(l 0 wt) 


(37) 


By knowing the character of the noise mixed with the measurement operator, the referee can change the j3s,t values 
and in this way lip) remains proportional to Tr(ITp). If he don’t know that, then the error will be proportional to 
Tr[(IT" - W)p\, where, 


/ / rri 11 rr\ rri rji 

Ts -Ts 

s,t 

and, Ps,t are calculated according to the equation [ 2 ]. 

With noisy input states and measurement operator 


W” -W = Y,Ps,t 


(38) 


lip) = liu, p,rrii,rn 2 ,ni,n 2 , p) = 7 ( 1 ,1,0,0,0, 0, p) 


+Tr 


P{W''- V IdsitTj 


s.t.a.b 


(39) 


So noise (both in preparation as well as measurement) can, in principle, degrade the quality of measurement device 
independent implementation of the entanglement witness operator. 
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B. Expression for ImodiPAB) 


Our original witnessing function was, 


I{pab) = 

S,t 


|$+) ($+1 0 |$+) ($+|)(ts 0 PAB ® Wt) 


and the modified form of that is, 


Iraod{pAB) = ^ /3s,t C Tr {PaoA ® PbBo){'^s ® PAB ® ^t) 

s,t 


where the modified measurement operators are given by, 

PaoA ® PbBo = MiM2COS^(Ai)cos^(A2) |$+) ($+1 0 |$+) ($■'■1 
+P1P2 cos^(Ai) cos(A 2) sin(A2) |<I)+) ($+| 0 ^ |' 1 '“) (<i)+| + |$+) ^ 

+P 1 P 2 cos2(Ai)sin2(A2) |$+) ($+| 0 l^-”) (^'"l 
+ m(l-M2) 1^ + ^ ($+1 0 I 

+P 1 P 2 cos(Ai) sin(Ai) cos^(A 2)^ |'1'“) ($+1 + |$’'‘) ^ 0 !<&+) (<i)+| 

+P 1 P 2 cos(Ai) sin(Ai) cos(A 2) sin(A2) 

( l^--) ($+1 + |$+) (^--l ) 0 ( l^--) ($+1 + |$+) (^--l ) 

+P 1 P 2 cos(Ai) sin(Ai) sin^(A 2 )^ ($■*■! + 1$+) (4'“| ^ 0 |\E'“) ('i'“| 

+ ^i(l ~ ^ 2 ) cos(Ai)sin(Ai)( 1^^) ($+1 + |$+) 0I 

+P 1 P 2 sin^(Ai) cos^(A 2 ) |\E'“) ('i'“| 0 1$+) ($+| 
+P 1 P 2 sin^(Ai) cos(A 2 ) sin(A 2 ) I'E'”) ('l'“| 0 ^ |'1'“) ($+1 + |$+) ^ 

+P 1 P 2 sin^(Ai) sin^(A2) (d/”! 0 [d/'”) (d/”! 
+ Ml(l^ (g)I[+ 11-^i^cos2(A2)l0 |$ + ) ($+1 

+-^^—cos(A 2) sin(A2)10 ( l^--) ($+1 + |$+) 1) 

+ -^^—^^sin2(A2)l0 1^-) (^-| + 

4 Id 


Therefore, from Equ. (40), (41) and (42) we have, 


fmod(p^B) — C 
■ C0S^(A2) 


„ 2 /'a ^„„„ 2 /'A ^ I /^ 2 ) 2 


P 1 P 2 COS^(Ai) C0S^(A2) + 


COS^(Ai) 


M2(l - Atl) _2/A ^ , ( 1 - Ml)(l - P 2 ) 

16 


I{pab) + additional term 


(40) 


(41) 


(42) 


(43) 


Thus, both multiplicative and additive noises are present in the system. 
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Now, using Equ. (27) and (43), we have [ for the choices /ii = /r 2 = /i and Ai = A 2 = A ], 


ImodiPAts) — Q E Ps,t [2a3^sin(2A)a;t + 2ai/r cos(2A)xt 


2a2/iyt — 2ai/i sin(2A)zi 


+2a3fi cos(2A)zt + 2bifi{cos{2A)xs + sin(2A)zs) — 2&3^(sin(2A)a;s — cos(2A)zs) 
-2b2pys - 2c33/x^ sm^{2A)xsXt + Ci3/x^ sin(4A)xsXt - Csi/x^ sin(4A)a::sa;t 
+2cii/r^ cos^(2A)a;sa;t - 2c23/r^ sin(2A)ysa;t + 2032//^ sin(2A)a;sy4 
—2c2ip,^ cos{2A)ysXt — 2ci2/x^ cos{2A)xsyt + 2ci3/x^ sin^(2A)zsa^t 
+2c3i/r^ sin^(2A)a;s-2:t + Cny^ sin(4A)2:sa:t + sin(4A)2:sa;t — Cuy^ sin(4A)a;sZt 

-C33y^ sin(4A)a;s-2;i + 2c3iy^ cos^{2A)zsXt + 2ci3/x^ cos^(2A)xsZt + 2c22P^ysyt 
- 2 ci 2 y^ sui{ 2 A)zsyt + 2c2i/i^ sin(2A)ys2;t - 2c32At^ cos{ 2 A)zsyt - 2 . 023 ^^ cos{2A)ysZt 
-2ciiy^ sin^(2A)zs2t + Ci3/x^ sin(4A)zs2t - C3i/x^ sin(4A)2;s2t + 2033^^ cos^(2A)zs2t + 2 


( 44 ) 





